Temperature changes when adiabatically ramping up an optical lattice 
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When atoms are loaded into an optical lattice, the process of gradually turning on the lattice is 
almost adiabatic. In this paper we investigate how the temperature changes when going from the 
gapless superfluid phase to the gapped Mott phase along isentropic lines. To do so we calculate 
the entropy in the single band Bose-Hubbard model for various densities, interaction strengths and 
temperatures in one and two dimensions for homogeneous and trapped systems. Our theory is 
able to reproduce the experimentally observed visibilities and therefore strongly supports the view 
that current experiments remain in the quantum regime for all considered lattice depths with low 
temperatures and minimal heating. 



PACS numbers: 



I. INTRODUCTION 



During recent years substantial progress has been 
made in cold atom experiments. Among the greatest 
achievements arc the experimental realization of the su- 
perfluid to Mott insulating transition with bosonic atoms 
in optical lattices [l| and the BEC to BCS crossover in 
fermionic quantum gases 0. Due to the good tunabil- 
ity the excitement to use cold atomic gases as a so-called 
quantum simulator has been hard to temper. Such quan- 
tum simulators are systems which mimic - in a control- 
lable and ultraclean way - simple strongly-interacting 
models such as the Hubbard model. With ultracold atom 
experiments the dream of using such quantum simula- 
tors to obtain new insight into the long-standing prob- 
lems of other research areas seems 'almost' feasible. One 
of the notorious challenges of condensed matter physics 
that might be solved in such quantum simulations is the 
fermionic Hubbard model [1, [j| , believed to be relevant 
for high temperature superconductivity 

Beside the use as quantum simulators, the unprece- 
dented tunability of cold atomic systems makes them 
very promising candidates for quantum computers. The 
experimental production of entanglement has been a 
large step forward into performing quantum computa- 
tions [1]. 

The prerequisite of applying the quantum simulator 
to unknown physics, is a complete understanding of the 
present experiments. One would expect full agreement 
between the current generation of bosonic experiments 
and theory, since the properties of the underlying ho- 
mogeneous Bose-Hubbard model are rather well studied. 
However, our understanding of the present experiments 
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is far from complete. The interpretation of the results is 
mainly complicated by two points: 

(i) The presence of an external trapping potential 
can cause the spatial coexistence of the super- 
fluid and Mott phases. This replaces the quantum 
phase transition by a gradual crossover, which can 
obscure the characteristic signal of the quantum 
phases. 

(ii) A change in temperature due to adiabatic or non- 
adiabatic origins can also hide the signature of the 
quantum phase transition, replacing it by a thermal 
transition. Considerable heating in current experi- 
ments would cast serious doubt 0, H, [§| on former 
interpretations. Knowledge of the temperature in 
an optical lattice is therefore highly desirable, but 
whereas the temperature of a weakly-interacting 
bosonic gas in a parabolic trap can be measured 
accurately, no reliable temperature measurement 
exists in the presence of a deep optical lattice. 

In this manuscript we reexamine the interpretation of 
recent experiments of bosons in optical lattices focusing 
on the effect of temperature in the presence of a trap- 
ping potential. We determine the lowest temperature 
that can be reached by loading the atoms adiabatically 
into the optical lattice. It is reasonable to assume that 
such processes are isolated and close to thermodynamic 
equilibrium at all times, since ramping up the lattice is 
typically slow (~ 16 ms per increase E'r in lattice laser 
intensity) compared to the tunneling rate of the bosons 
1 kHz, except for the deepest lattices) [l^. This as- 
sumption gives a lower bound for the temperature which 
can be achieved in the experiments without applying fur- 
ther cooling techniques. 

We compare the results of our approach with the 
measured experimental interference patterns over the 
whole parameter range. In contrast to previous work 



2 



0, H, [H, [3 1 ^6 '^^^ take the full range of realistic heights 
of the lattice potential and the full trapping potential into 
account and our calculations are approximation-free and 
based on first principles. 

After introducing the theoretical description of the 
quantum gas in section [ill discuss the Quantum 
Monte Carlo methods we developed in section Hill This 
section can be skipped by the non-expert. In section 
IIVI we present results for the entropy in one dimension. 
We compare homogeneous systems with systems in the 
presence of realistic trapping potentials across the super- 
fluid to Mott-insulating transition and we also study the 
Tonks-Giradeau gas. We make contact with analytic ap- 
proximations where possible, and we study the validity 
of the local density approximation. Section |Vl is ded- 
icated to the entropy in two-dimensional homogeneous 
and trapped systems. We study the influence of the trap- 
ping potential, density, temperature, and of the actual 
experimental procedure on the visibility in section IVII 
We also compare the visibility obtained computationally 
with experimental data. 



II. THEORETICAL DESCRIPTION 

In order to prepare ultracold atoms into an optical 
lattice, the atoms are confined in an external trapping 
potential and cooled till they Bose-Einstein condense. 
Thereafter, additional lasers forming an optical lattice 
are turned on. One gradually increases the intensity of 
these lasers, slowly enough to remain in a low-energy 
state but also fast enough such that external influences 
are small. 

Adiabaticity can be checked by reversing the process, 
i.e. decreasing the intensity and comparing the state at 
a certain lattice height shows the same experimental sig- 
nature as the corresponding state while switching on the 
lattice. Such experiments indicate that the loading of the 
atoms might to a good approximation be considered as 
adiabatic. 

Atoms loaded in a sufficiently deep optical lattice are 
described by the Bose-Hubbard model |14| . 



transition, but it becomes approximative for very low re- 
pulsion or high temperatures. 

The external trapping potential is included using = 
WcT^, where Vc describes the strength of the trapping 
and r the distance to the center of the lattice. If red 
detuned lasers are applied for creating the optical lat- 
tice potential the focusing of the lasers gives rise to an 
additional confinement. Under proper alignment of the 
two trapping potentials, the total trapping is given by 
e{i) = 771^^2/2 [H with 



Lo^ =Lo^ + 8Vo/{mw^), 



(2) 



where w is the waist of the lattice laser beam and Vq 
its intensity expressed in single-photon recoil energies, 
Eji = h^/2mX^. Here A is the wavelength of the lattice 
laser beam and m is the mass of the atoms. We took the 
waists isotropic and neglected corrections of the order of 
The second term in Eq. ^ dominates already for 
moderately strong lattices. 

The homogeneous Bose-Hubbard system shows a quan- 
tum phase transition from a superfluid to a Mott insu- 
lating state when the filling is commensura te [ l6l . The 
transition occurs at {U/t)c = 3.28(4) Il8l| in one- 
dimensional systems, and at {U/t)c = 16.74(1) [H, ^ 
in two-dimensional systems. Whereas in the superfluid 
state a continuous excitation spectrum exists at low en- 
ergies, the Mott-insulator is characterized by a gap just 
above its ground state. 

Analytically the entropy in the Bose-Hubbard model 
has been studied for non-interacting [13, [2l| and weakly 
interacting [l^.[l3| bosons. The strongly interacting limit 
for the Mott-insulator and the Tonks regime have been 
considered for homogeneous [H, [2l| and trapped systems 
[H, [111 . In three dimensions, the entropy deep in the su- 
perfluid and Mott insulating phases was calculated using 
effective masses |23| . Information on the trapped system 



was then obtained using the local density approximation 
(LDA). The Bose-Hubbard model can only be solved an- 
alytically in these limiting cases, where one is very deep 
in the superfluid or Mott insulating phase. Close to the 
phase transition numerical tools have to be employed, 
such as the QMC simulations performed here. 



H 



U 



^n,{n, - 1) + 



(1) 



The notation refers to the sum over nearest- 

neighbor sites only. Bosons are created on site j by the 
operator 6^ and the number of bosons on site j is counted 
by the number operator rij. The kinetic term describes 
hopping of the bosons with tunneling amplitude while 
the on-site repulsion has strength U . We will work in the 
canonical ensemble with a constant number of particles, 
A^. The linear system size is L and we work in dimen- 
sions d — 1,2. We restrict ourselves to a single-band 
Bose-Hubbard model, which is sufficient at low temper- 
atures and lattice intensities around the superfluid-Mott 



III. METHODS 

In this chapter we discuss a number of methods to 
determine the entropy 



S{(i) = (}{E -F) = I3E + ln(Z) 



(3) 



where E is the total (internal) energy of the system, 
f3 = 1/T the inverse temperature and F = —\n{Z)/(3 
is the free energy. We use /cb = 1. The main challenge 
is an accurate calculation of the partition function. It 
turns out that a combination of two methods discussed 
below gives the best results. Both are accurate, and do 
not involve any fitting nor noisy derivatives of numeri- 
cal data, but they are efficient in a different temperature 
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range. The results of both methods have been checked 
against each other for consistency. The fiat histogram 
methods (Sec. IIIISP works best at high temperatures, 
wliile the thermodynamic integration method (Sec.|TTTC| 
works better at low temperatures. This chapter is in- 
tended for the technically oriented readers and it is not 
necessary to read it in order to understand the discussion 
on the results in the next chapters. 



A. The canonical worm algorithm 

All our simulations employ a canonical worm algo- 
rithm. A worm algorithm [2J| in the path-integral rep- 
resentation is a quantum Monte Carlo algorithm where 
the decomposition of the partition function, 







Z = Tr^ / dtn I dtn-^i 



dti 



(4) 



p-tiHay^-{t2-tt)Ha . . . g-(^„-t„-l)ffo-^/g-(/9-^,^)ffo ^ 

is sampled indirectly by making local moves in the Green 
function sector, which is the extended configuration space 
of world lines with two open ends. Simulating the Bose- 
Hubbard model, we choose as diagonal part iJo the po- 
tential energy, while the hopping terms are the pertur- 
bation V. 

In a canonical worm algorithm the operators of the 
equal-time Green function hi{T)h^-{T) are propagated si- 
multaneously. The extended partition function we sam- 
ple reads 



Z, = Tr[r ((6,(r)6](r)+h.c.)exp(-/3H) 



- E E E 

«=0{|ji)},...,{|»„+i)} 
Wc{.)dti ■ ■ - dr ■ ■ ■ dtn, 

where the terms Wc(-) denote 



13 i-t 



(5) 



Hn\v\ 



-(t2-tl)B2 . . . 



. . „-(t„ + i-t„)B„ 



(in+i\V\ii)e' 



^(6) 



with Ei = (ii|i7o|*i)i and we have introduced sums over 
a complete basis set between any two off-diagonal opera- 
tors. The terms We are all positive and can thus be used 
as weights in a Monte Carlo sampling. 

An efficient updating scheme has been presented in 
Ref. H^, [11] and allows the straightforward computation 
of the kinetic and potential energy, density, compressibil- 
ity, equal time Green function, and superfluid density. 
However, the partition function is not a thermodynamic 
average and is harder to compute. 



B. Flat histogram methods 

The goal of a flat histogram method is to obtain a den- 
sity of states p{X) (where the coordinate X in classical 
simulations is usually the energy) directly by a random 
walk in A"-space instead of performing a canonical sim- 
ulation at fixed temperature. By sampling each value of 
X with a probability Q{X) oc 1/ p{X) we obtain a flat 
(constant) histogram H{X)p{X)Q{X) = const. 

In the Wang-Landau sampling scheme [l^], a crude 
guess for Q{X) is itcratively updated until it converges 
by a multiplicative factor /. During consecutive Wang- 
Landau iterations. / is reduced according to / 
when the current histogram H{X) is considered to be 
sufficiently flat. Then the histogram is reset, we have a 
more accurate estimator for the density of states, and the 
sampling restarts with the smaller /. The convergence of 
the scheme was proven by Zhou and Bhatt [1^ . In partic- 
ular, they showed that the minimum number of steps in 
each Wang-Landau iteration should scale as 1/ V/- The 
generalization of the Wang-Landau scheme to quantum 
systems was discussed in Ref. [l^. Here, we generalize 
Eq. (4) of Ref. to the path-integral formulation: 

Za = Trexp(-/?(ifo-AF) 

= TrV / dtn / dtn-i--- / dtie-^'"°XV 

g-(t2-tl)Ho . . . Q-{t„~t„-i)Ho )^Y^-{P~t„)Ho 



= E5(")^"- 



(7) 



71=0 



The expansion order n corresponds to the number of 
kinks (particle hoppings) present in the path integral 
representation of a configuration. The original partition 
function (which is a function of the inverse temperature 
(3) can be found back by setting A = 1. The density of 
states g{n) corresponds to 

E E r r^(-)^^i" 

(8) 

where W denotes the weight of a diagonal configuration, 
W{.) = e-*i^Mn|^l*2>e-(*^^*i)^^ ... 



-(i„-t„-i)-E„ 



in\V\ii)e 



(9) 



In such configurations all world- lines are continuous, and 
it occurs during the Monte Carlo run when the two open 
ends (worms) cancel each other on the same site and 
imaginary time. 

Using the canonical worm algorithm [H, [2^ , the den- 
sity of states can be obtained as follows : A single Monte 
Carlo step is defined from a diagonal to a new diagonal 
configuration and has an acceptance factor q' . Taking 
the density of states into account, the acceptance factor 
should be modified to 



q{x y) = min [1, g{x)q'/g{y)] 



(10) 
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where g{x) is the density of states corresponding to the 
expansion order of the old configuration x. When the 
expansion order of the new configuration is larger than 
a predefined maximum expansion order, the update is 
rejected. 

When the Wang-Landau iteration is finished, we can 
obtain the partition function for all values of A. For 
the Bosc-Hubbard model that means that we obtain 
a whole set of partition functions through the scaling 
l3t ptX,U/t U/Xt (thus PU is constant). A trap 
would also rescale as ^ ^ XV, which makes this scaling 
less useful in the trapped case and we use it just to obtain 
values at A = 1. If we had worked in the grand-canonical 
ensemble, the chemical potential would scale analogously, 
and we lose all control over the particle number. Here we 
see a distinct advantage of the canonical ensemble over 
the grand-canonical ensemble. 

The normalization of g{n) is fixed by calculating the 
partition function Zjv in the canonical ensemble for the 
case without hopping [t — 0). 



C. Thermodynamic integration 

The second method we discuss is the thermodynamic 
integration method. We choose a set of inverse tempera- 
tures 



where we sample all configurations a at the temperature 
f3j . In the canonical worm algorithm we have to measure 

. (l-f (15) 

with A/3 — (3j — (3j-i and Ed the time averaged potential 
energy of the configuration. We can thus compute the 
partition function at Pj if we know the partition function 
at f3j^i. 

The accuracy of the scheme depends on the overlap 
between the system at (3j and the one at Pj-i- If the 
overlap is small, the error on the partition function will 
increase rapidly and propagate systematically on to lower 
temperatures. In particular the first term Pi should be 
chosen sufficiently close to zero, since there are only con- 
tributions if the expansion order is zero. The fiuctuations 
in the (diagonal) energy in Eq. ([T5|) are exponentially 
hard to control. We therefore choose our set of values 
of P such that ApEd < 1. At large values of U or large 
particle numbers, more /3-points are required. When we 
are close to the ground state things get easier since there 
energy fluctuations are suppressed. In the limit that A/3 
is infinitely small, the scheme reduces to an energy in- 
tegration. Note however that the scheme remains exact 
when A/3 is finite. 



Po=0<Pi<P2<...<Pn. 

Then the partition function can be written as 



In Zi3^ = IhZq + In 



(11) 



(12) 



The partition function at infinitely hot temperatures 
Zo = ^/3o=o can be found by solving the combinatorial 
problem of placing N bosons on L lattice sites. 



In Zq ~ In 



L + A^-l 
N 



L+N~l 



N 



= E Inj-^lnj. (13) 



The ratios in Eq. ( [T^ can be estimated through the 
weights introduced in Eq. ( [5]) , 



(14) 



D. Numerical strategy 

Although we tried a number of alternatives, none of 
them were satisfactory. We briefiy make some remarks 
about them. We tried to use 



dS_ 

dU 



dT 



(16) 



and integrate the density fiuctuation with respect to in- 
verse temperature, but the scatter of the data was much 
bigger than the trend-line, which made this approach pro- 
hibitive without an adequate fitting procedure. Integrat- 
ing the specific heat 



S{p) 



1//3 



cv{T')/T'dT', 



(17) 



after differentiating the fitted curve through the en- 
ergy was used in Ref. [H, [SO, [IH. However, the divi- 
sion by the temperature is misbehaving at low T and 
the specific heat computed via the fluctuation formula 
cv = /3^((£'^) — {E)"^) is a quantity that converges 
slowly in the Monte Carlo simulation. We have also 
combined a grand-canonical directed loop algorithm in 
the path-integral representation (s^ ) with a quantum 
Wang-Landau rcweighting scheme, though the fact that 
the Wang-Landau reweighting also changes the density, 
made this approach very cumbersome. A better attempt 
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TABLE I: Comparison between the flat histogram ('QWL') 
and direct integration ('chain') method for a homogeneous 
one- dimensional Bose-Hubbard system with f//f = 2 and A'' = 
45. 



pt 


Sqwl 


5'chain 


0.5 


44.2(2) 


44.08(3) 


1.0 


28.1(2) 


27.83(3) 


1.5 


19.0(3) 


18.6(1) 


2.0 


13.6(4) 


13.2(1) 


3.0 


7.9(4) 


7.6(6) 



TABLE II: Comparison between the flat histogram ('QWL') 
and direct integration ('chain') method for a trapped two- 
dimensional Bose-Hubbard system with U/t — 100 and Vc/t — 
2.5. 



I3t 


Sqwl 


5'chain 


0.2 


36.9(3) 


36.7(2) 


0.5 


12.2(3) 


12.1(2) 


0.7 


7.2(4) 


7.4(4) 


1.5 


1.9(2) 


2.0(4) 



w 0.5 
0.4 




5 10 15 20 

T/t 



FIG. 1: Entropy per particle as a function of temperature for 
a homogeneous one-dimensional Bose-Hubbard system with 
U/t = 100, A = 40 and L = 50 (n = 0.8). The plateau hints 
at the existence of a gap in the excitation spectrum. 



IV. ISENTROPIC LINES IN ONE 
DIMENSIONAL SYSTEMS 



was developing a canonical directed-loop algorithm in the 
stochastic series expansion [ssj and combining it with a 
quantum Wang-Landau reweighting scheme. This ap- 
proach has the advantage that one obtains all tempera- 
tures down to the one corresponding to the pre-chosen 
cut-off length at once. The drawback is that in a SSE 
representation the large values of U are also sampled, 
requiring extremely large orders even for moderate tem- 
peratures. 

We obtained the highest accuracy by combining 
the two methods outlined in detail above. For high 
temperatures, up to /3t « 0.5, the combination of 
the canonical worm algorithm with the flat histogram 
(QWL) scheme is best and fast. For larger (3t, and 
since we arc interested in the entropy for a very large 
number of pt over a relatively small temperature range, 
thermodynamic integration is best. This method should 
only be used close to the ground state. Otherwise the 
fluctuations in Eq. are hardly controllable which 

might lead to large systematic errors. 

We compare the accuracy of both methods for a 
homogeneous one-dimensional system and for a trapped 
two-dimensional system in Table U and Table |TT1 re- 
spectively. The data and the error bars shown in 
the tables do not reflect the full computational cost, 
since for the propagation method we need many more 
intermediate values of /3 which are not shown in the 
tables. Similarly, for high (3 the flat histogram method 
requires only a single Monte Carlo run, but the cost 
scales exponentially with /3 (and the mmrber of particles). 



We will start our discussion of the one-dimensional case 
with a homogeneous lattice and compare our results with 
analytic results in limiting cases. We will then gradu- 
ally make the discussion more realistic (and more com- 
plicated) by including the parabolic confinement. The 
first step is discussing the entropy in a system of constant 
quadratic trapping, Vc/t = const.. The second step is the 
case of an external parabolic trapping potential which 
is further strengthened by the focus of the lattice laser 
beams as is currently done in most experiments. We will 
compute the entropy with these two confining potentials 
starting in the superfluid and going to the Mott-insulator 
or the Tonks-Giradeau gas. We will also check the qual- 
ity of a numerically based local density approximation 
and will find that in its regime of validity the speed-up 
in computing the entropy of trapped systems is consid- 
erable. 



A. Homogeneous system 

The dependence of the entropy on the temperature 
is closely related to the energy spectrum. Weakly in- 
teracting superfluid systems have an continuous energy 
spectrum. Therefore the entropy rises continuously 
with temperature. Mott insulating states have an 
energy gap just above the ground state and entropy 
is exponentially suppressed up to temperatures of the 
order of ksT/U « 0.1 [ll,[23. Above the gap the bands 
of excited states lead to a finite entropy if temperature 
is high enough. Strongly interacting incommensurate 
systems also have gaps in the spectrum, but at sub- 
stantially higher energies, signaled by a plateau in the 




0.5 1 1.5 2 2.5 3 

n 



FIG. 2: Entropy per site as a function of the filling factor for a 
homogeneous one-dimensional Bose-Hubbard system. QMC 
data are shown by the symbols for U/t — 12 and U/t = 
2 when jSt = 2 and L = 50. The lines are guide to the 
eyes. Comparison is made with the atomic approximation 
at different temperatures where only particle - hole (Ip-lh) 
contributions are taken into account. 



entropy in Fig. [TJ 

The dependence of the entropy per site on the filling 
is shown in Fig. [2] for different interaction strengths at 
a moderately low temperature (3t = 2. For weakly in- 
teracting systems {U/t = 2) the entropy reaches a max- 
imum at half filling and decays monotonously for higher 
filling factors. Mott regions in strongly interacting sys- 
tems {U/t = 12) appear as strong dips in Fig. [21 where 
the entropy is exponentially suppressed because the tem- 
perature is well below the gap. We also make a compari- 
son with the atomic limit approximation, where only sin- 
gle particle-hole excitations are taken into account. The 
agreement with the Monte Carlo data is only qualitative, 
because the atomic limit misses higher order particle- 
hole excitations which are important for U/t ~ 12. The 
atomic limit is expected to work well near commensura- 
bility, but already for densities 0.9 and 1.1 the deviation 
with Monte Carlo is considerable. The atomic limit ap- 
proximation also incorrectly predicts a symmetric curve 
around n = 1. 

Figure [3] shows the entropy per site for a system of 40 
particles on a lattice of 50 sites as a function of the inter- 
action strength U/t and inverse temperature t/T = (it. 
Since we are away from integer filling we are always in 
the superfluid phase. At intermediate and high temper- 
atures /3t < 3, the entropy decreases when increasing the 
interaction strength. In contrast at lower temperature, 
P't ^ 3, we see a non-monotonic behavior of the entropy 
with a minimum close to the critical interaction strength 
in a commensurate system. The same qualitative behav- 
ior was observed for a smaller lattice of L = 20 sites. 
The reason is the presence of the nearby quantum phase 
transition to a commensurate Mott state. The mass of 
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1 2 3 4 5 



t/T 

FIG. 3: Entropy per site of a homogeneous one-dimensional 
Bose-Hubbard model. There are 40 particles and 50 sites at 
incommensurate filling. The system is always in the superfluid 
regime. The error on the entropy is less than two percent in all 
cases. This representation of the entropy allows to read off the 
final temperature easily when adiabatically ramping up the 
lattice. For instance, if the initial temperature is /3£ = 2 for 
U/t = 2 then the final temperature is pt = 1.3 for U/t = 10. 

the quasi-hole decreases when going away from the tip 
of the Mott lobe at {U/t)c to higher values of U/t. The 
quasi-hole absorbs more entropy and the entropy thus 
increases. Moving along isentropic lines we observe in 
units of t some mild heating when the initial tempera- 
ture is reasonably high, /3t < 3, but cooling if the initial 
temperature is sufficiently low, (3t > 4. 

The commensurate case in Fig. 2] shows the same be- 
havior as Fig. [3] for low values olU/t < {U/t)c when we 
remain in the superfluid regime. However, when the Mott 
state develops, i.e. U/t > {U/t)c, the gap in the spectrum 
opens and the entropy at constant temperature decreases 
considerably. Along adiabatic trajectories the tempera- 
ture in units of t shoots up near the transition point. 
This has been discussed for the homogeneous system in 
Ref. [2l[, from which the authors deduced that the tem- 
perature in present experiments must be of the order of 
U. However, Rey et al. pointed out that in the pres- 
ence of a parabolic confining potential less heating occurs 
in hard-core bosonic systems than in a homogeneous sys- 
tem. The next section addresses the same question for 
soft-core bosons. 



B. Entropy distribution in a constant parabolic 
trapping potential 

In the presence of an external trapping potential spa- 
tially separated quantum phases can coexist. In Fig. [5] 
we show two density profiles in the presence of a trapping 
potential. The first is for a superfluid state (Fig. [5^)) in 
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FIG. 4: Entropy per site of a homogeneous one-dimensional 
Bose-Hubbard model at commensurate filling. There are 50 
particles and 50 sites. For U /t < 3.28 [l7l.[l8l| the system is in 
the superfluid phase in the thermodynamic limit at T = 0, for 
larger U/t a, gapped Mott state is formed. The error on the 
entropy is of the order of 0.05, making it impossible to observe 
the exponential decay of the entropy in the Mott state. 



which the density closely follows the form of the trap- 
ping potential. The large local variance k = (nf) — (rii)'^ 
demonstrates the number fluctuations in the superfluid 
state. The second density profile in Fig. [5Jd) is for a state 
in which the central inhomogeneous region is surrounded 
by a Mott-insulating shell with commensurate filling and 
an outer incommensurate region. The variance shows a 
clear suppression of the density fluctuations in the Mott- 
insulating regions. 

The coexistence of spatially separated quantum phases 
can be understood in terms of a site-dependent chemi- 
cal potential, the so-called local density approximation 
(LDA). Physical quantities are determined by using on 
each site the results obtained for a homogeneous system 
with the corresponding chemical potential. The site- 
dependent effective chemical potential provides a scan 
through the phase diagram. This approximation has 
been shown to work nicely for such quantities as the 
density or the variance in regions where the trapping po- 
tential varies slowly [s^, [sH]. However, the LDA breaks 
down for steep trapping potentials and near the edges of 
Mott plateaus where numerical simulations are necessary 
to obtain reliable values [s^l . To get a better understand- 
ing of the entropy distribution in an inhomogeneous sys- 
tem we developed a canonical and improved variant of 
the LDA, dubbed iLDA: we first calculate the exact den- 
sity profile using a full QMC simulation for the trapped 
simulation. Then we take for every site the entropy from 
a homogeneous run corresponding to that density. This 
variant has the advantage that we start from the exact 
density profile, taking into account the rounding near the 
edges of the Mott plateaus due to the finite gradient of 
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FIG. 5: Density profile (red line, top curve in the center), 
density variance (green line, middle curve in the center) and 
entropy profile (blue line, lowest curve in the center) for a 
one-dimensional trapped system for different parameters, a) 
U/t = 2, inverse temperature I3t — 2, system size L — 128, the 
trapping potential Vc/t = 0.00829, and the density = 80. 
The total entropy calculated by the LDA approximation is 
5lda = 10.75 and by QMC simulations Sqmc = 9.8(1). b) 
U/t = 12, inverse temperature l3t = 2, system size L = 128, 
the trapping potential Vc/t = 0.00829, and A'^ = 80 particles. 
Sqmc = 9.2(1) from QMC and Slda = 8.3 from LDA. 



TABLE III: Entropy comparison between iLDA and exact 
Monte Carlo results for different one-dimensional system pa- 
rameters. The trapping potential for the parameters in the 
third row is Vc/t = 0.00829. 



parameters 


S by QMC 


S by iLDA 


Fig. [5l left panel 




9.8 ±0.1 


10.8 


Fig. [5l right panel 




9.2 ±0.1 


8.3 


U/t = 6, f3t = 1,N = 


80 


22.7 ±0.1 


21.8 


Fig. [7], U/t ^12 




7.9 ±0.2 


6.6 



the trap. We tested the approach by comparing the total 
entropy of the trapped system calculated as the sum of 
the single sites to full numerical simulations. We found, 
as shown in Tab lllll that the iLDA can capture the qual- 
itative trend of full calculations, but cannot reproduce 
the exact values. 

Using iLDA we obtained the entropy profiles for the 
density profiles of Fig. [S] In a system where superfluid 
and Mott-insulating regions coexist, we clearly see that 
the Mott-like regions are not able to accommodate en- 
tropy and the whole entropy is in the superfluid regions. 
If the whole system is superfluid (Fig. [5] (a)) the entropy 
varies only slowly from site to site and shows maxima at 
the flUing close to n « 1/2. Since the flUing is larger than 
one in the center a dip in the entropy profile develops. 

In Fig. [HI we show the entropy calculated in a full QMC 
calculation for a constant trapping potential. In Table HVl 
we show the values of temperatures following an isen- 
tropic line as extracted from the data. The concept of 
'adiabatic heating' is complicated by the different energy 
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t/T 

FIG. 6: Entropy per particle of a one-dimensional Bose- 
Hubbard model with constant trapping Vc/t = 0.00829. 
There are 80 particles and 128 sites. The magnitude of the 
errors is approximately 0.1 — 0.2. The roughness of the lowest 
isentropic line is within the error bars. 
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FIG. 7: (a) Density profiles and (b) variance along the isen- 
tropic line of S = 20 in Fig. [O] See Table IIVI for the corre- 
sponding temperatures. 



TABLE IV: Values of U/t, pt and Uf) along an adiabatic line 
5 = 25 for the same parameters as in Fig. [S] See Fig. [7] for 
the corresponding density profiles. 



U/t 


(3t 


pu 


4 


1.10(2) 


4.4(1) 


9 


1.10(2) 


9.9(2) 


11 


1.00(5) 


11.0(5) 


13 


0.93(3) 


12.1(4) 


15 


0.80(5) 


12.0(8) 


18 


0.65(3) 


11.7 (5) 



TABLE V: Values of U/t, pt and PU along an isentropic line 
S = 13 for the same parameters as in Fig. |8] 



Vo[En] 


U/t 


Vc/t 


pt 


pu 


4 


6.13 


0.008 


1.62(4) 


9.9(2) 


5 


8.62 


0.011 


1.35(4) 


11.6(4) 


6 


11.90 


0.015 


1.10(4) 


13.1(4) 


8 


21.74 


0.028 


0.70(4) 


15.2(8) 


9 


28.57 


0.037 


0.61(2) 


17.5(6) 


10 


37.04 


0.049 


0.49(2) 


18.2(7) 



scales (and units) used in the literature. It is important 
to define with respect to which energy scale the temper- 
ature is measured as illustrated in Table [TVl We see that 
along isentropic lines in the superfluid phase {U/t < 10) 
temperature in units of the hopping t remains roughly 
constant. At higher temperatures there is some small 
heating in units of t, while at low temperatures we ob- 
serve a little cooling. 

This behavior can qualitatively be understood by look- 
ing at the density profiles and its variance along isen- 
tropic lines for different values of U/t, as shown in Fig.[7l 
Prior to the formation of a wide commensurate region, 
the temperature remains almost constant, for instance 
pt = 1.1 for U/t = 4 and U/t = 9 and /3t = 1.0 for 
U/t = 11. Only when a considerable volume percentage 
of the system turns insulating the incommensurate edges 
cannot accommodate the entropy anymore which results 
in a rise of the temperature, i.e. f3t = 0.8 for U/t = 15 
to pt = 0.7 for U/t = 20. We checked that this effect is 
seen for a wide range of initial temperatures. 

In units of the interaction strength U a cooling takes 
place. In particular, we see that for the chosen ini- 
tial temperature the final temperature with respect to 
U stays below the temperature for which excitations in 
the Mott insulator are created in a homogeneous system. 
Therefore, the Mott insulator is stable up to the consid- 
ered lattice height and the superfiuid regions take most 
of the excitations. 



C. Entropy distribution in a realistic parabolic 
trapping potential 

Figure [8] shows the entropy when the finite waist of 
the optical lattice beam is taken into account (and v^/t 
is no longer a constant but a function of the lattice laser 



9 




0.5 1 1.5 
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FIG. 8: Entropy per particle of a one-dimensional Bose- 
Hubbard model with varying trapping potential, taking 
the waist of the laser into account according to Eq. ([2]) 
with a magnetic trapping frequency ojo = Stt x 30 Hz 
and a laser waist w = 160/im for *^Rb atoms. Data 
points were calculated for Vb = 1, 2, . . . , IO-Er or for 
U/t = 1.9,2.9,4.3,6.1,8.6,11.9,16.1,21.7,28.6, and 37.0 and 
for A/3t ~ 0.1 or 0.2. There are 60 particles and 128 sites. 
The magnitude of the errors is a few percent. 



TABLE VI: Values for U/t, fSt, and PU along an isentropic 
line 5* ~ 20 for the parameters as in Fig. |9l See Fig. [lO] for 
the corresponding density profiles. 



Vo[Er] 


U/t 


tp 


Ul3 


1 


1.92 


1.70(2) 


3.2(1) 


4 


6.13 


1.15(5) 


7.0(3) 


5 


8.62 


1.00(3) 


8.6(3) 


6 


11.90 


0.85(5) 


10.2(6) 


7 


16.13 


0.70(2) 


11.3(3) 


8 


21.74 


0.55(4) 


12.0(9) 


9 


28.57 


0.46(2) 


13.1(6) 



intensity) . Table |V] gives the parameters for the tem- 
perature extracted along the isentropic Une = 13. We 
see that along this line the temperature in units of t in- 
creases while in units of U the temperature decreases. 
This means that the formation of the Mott-insulator is 
still possible starting at a low enough temperature, since 
the incommensurate regions take a lot of entropy. 

Figure [5] shows the entropy as a function of the temper- 
ature for different strengths of the optical lattice poten- 
tial. Compared to Fig.[5]the number of atoms is increased 
to iV = 140. This has the consequences that even at low 
temperature only small Mott-insulating regions can form 
and an incommensurate region survives in the center of 
the trap (Fig.[Tni). In Table IVll we show the temperature 
along the isentropic line = 20. In units of t we now 
get very strong heating, even in the superfluid phase, due 
to the high occupation which does not accommodate as 
much entropy as the low occupation region (cf. Fig[2|). 
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FIG. 9: Entropy per particle as a function of the inverse tem- 
perature /3 for difi'erent intensities Vb (expressed in recoil en- 
ergies) of the lattice laser. The system has N = 140 parti- 
cles on a lattice of L = 128 sites. We took a laser waist of 
w — 160^m. An isentropic line of 5 = 10 goes through 13 — 3 
for Vo = 1 and through /3 = 1 for Vb = 8. See Table for 
temperatures along an adiabatic line S = 20. 
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FIG. 10: (a) Density profiles and (b) density variance profiles 
along the isentropic line of S = 20 in Fig. |9] See Table IVll for 
the corresponding temperatures. 
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TABLE VII: Values for U/t, /St, and /SU along an isentropic 
line 5* ~ 50 for the parameters as in Fig. [9] 



Vo[Er] 


U/t 


pt 


pu 


1 


1.916 


0.80(2) 


1.5(1) 


3 


4.23 


0.60(2) 


2.5(1) 


4 


6.14 


0.50(2) 


3.1(1) 


7 


16.13 


0.27(2) 


4.4(5) 


8 


21.74 


0.22(2) 


4.8(4) 




12 3 4 

t/T 

FIG. 11: Entropy per particle of a one-dimensional Bose- 
Hubbard model with constant trapping Vc/t — 0.00829, ap- 
proaching the Tonks regime for large values of U/t. There are 
N = 15 particles in the system of size L = 50; the system is 
so dilute such that no Mott region is formed. 

However, in units of U we again find a temperature de- 
erease. 

One sees that the Mott-insulator is stable against the 
temperature change, since the temperature stays below 
O.IU. Since the Mott-insulating regions are small and 
a broad incommensurate region survives in the center, 
the Mott transition does not play a central role in the 
behavior of the entropy curves. 

Starting at a higher initial temperature (3t = 0.8 and 
following the isentropic line 5* ~ 50, the same qualita- 
tive effect of temperature increase in the units of t and 
decrease in the units of U can be seen in Tab. IVIII 
However, at U/t w 20, the temperature is still so high 
{ksT > Q.IU) that no clean Mott-insulating region can 
be formed. 



D. Tonks gas : one-dimensional trapped case 

When the potential energy between the atoms of a one- 
dimensional Bose-Hubbard model increases, the particles 
behave more and more like hard-core bosons. The limit of 
infinite repulsion and no multiple occupancies of a site, is 
called the Tonks- Giradeau gas. Quantities such as the en- 
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FIG. 12: Entropy per particle of a one-dimensional Bose- 
Hubbard model with changed trapping, approaching the 
Tonks regime for large values of U/t. There are A'^ = 15 
particles in the system of size i = 50 



ergy, average density, variance of the density can be com- 
puted accurately by assuming non-interacting fermions. 
Other quantities such as the density matrix map to a non- 
interacting fcrmionic density matrix up to a phase factor 
coming from the Jordan- Wigncr transformation. The ex- 
perimental detection of the Tonks gas has demonstrated 
one of the fundamental concepts of quantum mechanics, 
namely the absence of a clear meaning of statistics in 
one-dimensional systems. 

The Tonks regime has been observed with [l^] and 
without a lattice [s^. In the experiment with a lat- 
tice, the data was analyzed using fermionization, and 
good agreement with experiment was found in the region 
U/t > 5. The fermionization results were obtained at dif- 
ferent temperatures along adiabatic lines. Temperature 
rose from (3t = 2 for a lattice depth of Vb = 4.6£'r till 
fit = 0.77 for Vb = 12£'p{ as a consequence of the change 
in Vc/t when ramping up the lattice. Consistently, it was 
argued in Ref. 0] that the temperature in the Tonks gas 
in a trap was of the order of the hopping t. 

The Tonks problem was also studied by Monte Carlo 
simulations [33, [sll at a low but constant temperature 
fit ~ \. which is of the same order as the one used in the 
fermionization approach. The authors compared hard- 
core and soft-core bosons for a homogeneous lattice and 
for constant trapping. They found a gradual cross-over 
and found that the presence of a trap did not qualita- 
tively change the Tonks onset. 

The contrast in the experimental interference pattern 
was almost completely gone for the deepest lattices [2^ . 
Although this is consistent with a strongly repulsive su- 
perfluid gas, one might fear that similar patterns are pro- 
duced by either a Mott state or a thermal state due to 
a combination of soft-core bosons and an increased trap- 
ping depth. 
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Our analysis is again carried out in two steps. First, 
we calculated the entropy of a one-dimensional super- 
fluid in the very low density limit trapped in a constant 
parabolic trap. From Fig. [TT] we deduce that temper- 
ature remains remarkably constant. Thus, for a super- 
fluid in the low density regime, adiabatic processes are 
(almost) isothermal when the external trapping is con- 
stant and weak. Second, we make the simulation of 
the experiment more realistic. We numerically evaluate 
the Bose- Hubbard parameters using the tight-binding ap- 
proximation, and obtain the same parameters U and t as 
in Ref . [sj ■ In contrast to Ref . [sj , we now also calculate 
the trapping parameter Vc/t from the total axial trapping 
Wax = 27r X 60 Hz for all optical lattice depths. We 
find in Fig. [T^] some heating, even in the low-density su- 
perfluid phase. Along an adiabatic line similar to the one 
taken in Ref. [1^ , we find that the temperature increases 
from f3t = 2 at Vo = 5Er (U/t = 6.9) to /3t = 0.69 
at Vo ^ 12Eb, [U/t = 49.5). Thus, the temperature in- 
crease compares v ery well to the one calculated assuming 
hard-core bosons [2^ . 

Summarizing, we confirm that the experiment has in- 
deed reached the gradual cross-over toward the Tonks 
regime. The temperature remains of the order of the hop- 
ping i, even though a temperature increase of a factor of 
3 is foimd due to the change in confinement strength Vc/t 
when ramping up the lattice. It is exactly this increase 
in temperature that prevents the Mott domains from de- 
veloping, since for Vq = 12E^ {U/t = 49.5, Vc/t = 0.073 
) a broad Mott domain appears in the center of the trap 
for 15 particles and (it = 2. At a temperature jSt = 0.69 
the central density is 0.9 however. 



V. RESULTS IN TWO DIMENSIONS 
A. Homogeneous two-dimensional superfluid 

We again start our analysis in two dimensions with the 
homogeneous case. The phase transition from the super- 
fluid to the Mott phase occurs for a, n = 1 commensurate 
system at {U/t)c = 16.74(1) [ll,[23|. We calculate the en- 
tropy for a superfluid system close to commensurability, 
n = 0.8. In Fig.[T3]the dependence of the entropy on the 
temperature and the interaction strength is shown. Its 
main behavior is similar to the one in a one-dimensional 
homogeneous case as reported in Fig. [1] At fixed inter- 
action strength the entropy shows only a small increase 
for low temperatures. However, at a certain temperature 
it starts to increase strongly before it bends down again 
and a plateau is formed. This can as in the ID case be 
related to the underlying band-structure in which first 
excitations in the lowest band can be created. Above a 
certain temperature the corresponding gap in the energy 
band-structure causes an intermediate saturation before 
at even higher temperatures further bands can be excited. 

As a function of the interaction strength, the entropy 
shows - for constant but low temperature - a minimum 
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FIG. 13: Entropy per particle for a system of size L = 20 x 
20 in the superfluid phase, A'^ = 320. For a commensurate 
system, the transition happens at {U/t)c ~ 16.74(1) [TqI. [20l|. 
For low temperatures, we see an initial heating with increasing 
U/t, but around the transition point the presence of the Mott 
phase is felt and the system starts to cool, thanks to the 
lower effective mass. For larger values of U /t we see a further 
cooling for low temperatures since the Mott phase is far away 
and we go deeper inside the superfluid phase. The inset shows 
the entropy for low temperatures (same axes and symbols as 
in the main flgure). 
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FIG. 14: Double occupancies for the same system as in 
Fig. [I3l The curves bend down for low temperatures at low 
values of U/t. For large values of U/t a, minimum is reached 
around /3t ~ 0.5 while for larger values of j3t the curves bend 
slightly upwards with temperature. 



close to the superfluid to Mott-insulator transition point 
of a commensurate system (inset of Fig. [T5|) . As in the 
one-dimensional case this can be attributed to the effec- 
tive mass change of the quasi-hole which has its maxi- 
mum close to the phase transition point for the homoge- 
neous system. 

The density fluctuation shown in Fig. [T3] is consistent 
with the behavior of the entropy shown in Fig. 1131 us- 
ing the relationship of Eq. For inflnitely hot tcm- 
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TABLE VIII: Bose-Hubbard parameters chosen in the two- 
dimensional trapped system. 



U/t 


Vc/t 


5 


0.2 


10 


0.33 


15 


0.47 


20 


0.61 


25 


0.74 


50 


1.11 


100 


2.50 



peratures the density fluctuation is independent of U/t 
(not shown). For low values of U/t density fluctuation 
goes down monotonously with (3t. The normal-superfluid 
transition happens around pt w 0.30(5) in our system of 
density 0.8 for a small lattice of 20 x 20 and U/t^'i. This 
transition belongs to the Kosterlitz-Thouless university 
class, and was studied in detail for the commensurate 
case in Ref. 

For large values of U /t we enter the quantum (ther- 
mal) critical regime determined by the quantum critical 
point U/t = 16.74 [l9l.[20| and a minimum in the density 
fluctuations is reached. This is clearly seen for U/t = 15 
around f3t = 0.5. For larger values of U/t the minimum 
is reached at lower values of f3. After this minimum is 
reached, the density fluctuations go slightly up with (3t. 
The increase in the density fluctuations can be under- 
stood from the tendency of a dilute gas of vacancies (with 
respect to the n = 1 Mott state as a vacuum) trying to 
condense. In our canonical simulation, we will observe a 
tendency to increase the number of vacancies which will 
enhance pair formation. Thus theoretically the density 
fluctuations contain a lot of information about the sys- 
tem, but from the practical point of view the almost flat 
slopes in the quantum regime make this quantity a bad 
candidate for thermometry. 



B. The superfluid to Mott insulating transition in 
a parabolic trap 

In two dimensions we only consider trapping potentials 
taking into account the influence of the flnite waist of 
the lattice laser on our sample. We use the parameters 
shown in Table. Eml There are N = 200 bosons and the 
total system size is i = 20 x 20. These parameters are 
a compromise between increasing the trapping frequency 
while still obtaining meaningful results on a lattice of 
size L = 20 X 20. The fiUing in the system with TV = 200 
atoms for weak interactions is chosen such that it is close 
to n = 2 in the center for large U/t. 

The entropy dependence on temperature is shown in 
Fig. [in] for different values of U/ 1. For low values of the 
temperature the slope for the curves is very flat. For 
higher temperature a clear increase in the entropy in the 
system can be seen. For low temperature the curves for 
different U /t almost coincide and only the curve for very 
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FIG. 15: Entropy per particle as a function of inverse tem- 
perature for a set of different values of U /t with a trapping 
potential according to Table IVIIII The curves for the inter- 
mediate curves coalesce within error bars for fit > 1.5. The 
curve for U/t — 100 is significantly below the other curves 
(see however text) up to temperatures of the order jSt = 2. 
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FIG. 16: One-dimensional cut of the density profiles for y = 
L/2 depending on the x coordinate (labeled with 'i') along 
an adiabatic line S = 10. Error bars are smaller than the 
point size. See Table IVIIII for the corresponding values of 
the trapping potential and Table IIXI for the corresponding 
temperatures. 

strong interactions, i.e. U/t = 100, shows a considerably 
lower value of the entropy. 

The behavior of the entropy can be explained by con- 
sidering the density proflles of the system, shown in 
Fig. [in] along an adiabatic line S = 10. By increasing 
the optical lattice potential the central density goes down 
and we see the appearance of the n = 1 Mott-insulating 
region at U/t « 25. These regions can be identified in 
the variance proflle in Fig. [IT] For stronger interactions, 
U/t ^ 50, small Mott-insulating regions with n = 1 and 
n = 2 exist. For U/t = 70 the Mott-insulating regions 
cover already a large volume fraction of the system as 
clearly signaled in the variance proflle. The formation of 
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FIG. 17: Variance of the density along the x direction (la- 
beled with 'i') for y = L/2 along an adiabatic line S = 10, 
corresponding to the density profiles shown in Fig. 1161 Error 
bars are smaller than the point-size. See Table IVIIII for the 
corresponding values of the trapping potential and Table HXl 
for the corresponding temperatures. 



TABLE IX: Temperature along three adiabatic lines for the 
two-dimensional trapped Bose-Hubbard model with parame- 
ters according to Table IVlIIl 



U/t 


/3s=iot 






l3s=4oU 






5 


1.40(5) 


7 


0.70(2) 


3.5 


0.13(1) 


0.65 


10 


1.20(5) 


12 


0.60(2) 


6 


0.08(1) 


0.8 


25 


1.30(5) 


32 


0.50(2) 


12.5 


0.06(1) 


1.5 


40 


1.25(5) 


50 


0.40(2) 


16 


0.06(1) 


2.4 


50 


1.20(5) 


60 


0.35(5) 


17.5 


0.06(1) 


3 


70 


1.18(2) 


83 


0.25(1) 


17.5 


0.03(1) 


2.1 


100 


0.58(3) 


58 


0.18(2) 


18 


0.01(1) 


1 



the large Mott-insulating region causes the value of the 
entropy for the curve at U/t ^ 100 to he below the others 
at low temperature. 

The temperature evolution along isentropic hues for 
different initial temperatures is shown more quantita- 
tively in Table IIXI When starting from a low tem- 
perature, we see almost no heating in units of t up to 
U/t = 70. The initially low entropy can be distributed 
over the remaining superfluid regions. Measuring the 
temperature in units of U leads to a cooling of the system 
below the value oiU/t = 70. 

In contrast, for interaction strengths larger than U /t = 
70 almost the whole system is occupied by a commensu- 
rate region and the incommensurate regions have a neg- 
ligible volume fraction. The energy cost to generate an 
excitation in this situation corresponds in the bulk to 
the large interaction energy and at the boundaries to the 
large potential energy cost resulting from the steep trap- 
ping potential. Hence the entropy cannot be well accom- 
modated in the system and the temperature in units of t 
increases stronger than before when going to U/t = 100. 
In units of U the temperature first drops before it in- 
creases again or saturates for large U/t lattice potential. 



Starting with a low temperature ensures that the tem- 
perature remains low enough for the presence of the 
Mott-insulating regions. In contrast when starting from 
a hot temperature (the right-most example in Table IIXI 
is already in the normal state), there is heating in units 
of t and only weak cooling in units of U . The quantum 
degeneracy regime is never reached. Note that the de- 
scription by the one-band Hubbard model breaks down 
at such high temperatures. 



VI. INTERPRETATION OF EXPERIMENTAL 
RESULTS 

One of the standard experimental observation tech- 
niques consists of suddenly switching off the confinement 
and taking absorption images of the freely expanding gas 
after a finite flight time. The hereby resulting interfer- 
ence pattern is a reflection of the initial momentum dis- 
tribution, 
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where the factor w(k) is the Fourier transform of the 
Wannier function |14l | . The use of the momentum profiles 
(and also of the visibilities) has been subject to debate, 
since in this quantity different effects as heating or the 
loss of coherence by stronger interactions have the same 
consequence. Furthermore, it is very difficult to extract 
information on the superfiuid-Mott insulator transition 
point from these measurements in a trapping potential 
[ITI . [3^ . [sot - This is due to the spatially coexisting re- 
gions of superfluid and Mott-insulating character. A local 
measurement has to be implemented to obtain detailed 
information about the system. A first step into this di- 
rection has been taken by Foiling et al. [40| who mea- 
sured the density in thin slices. Experimental progress 
was reported for a local measurement of the density in 
Ref. [41|, and another proposal was made in Ref. [43 |. 
Such measurement techniques should be preferred over 
the visibility, which is only well suited for identifying the 
Mott and superfiuid phase far away from the transition 
point. In Figs. [TBI and [T7I we show the typical evolution 
of the density profiles (with Mott plateaus for strongly 
repulsive systems) and variance profiles along isentropic 
lines. 



A. Dependence of visibility on temperature and 
trapping potential 

Before we compare our results to the experimentally 
extracted quantities we would like to point out some fea- 
tures of the momentum distribution in a trapped systems 
at finite temperature. To do so we define the visibility V 
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FIG. 18: Visibility for a two-dimensional system L = 20 x 20 
as a function of U/zt at a high temperature of pt = 0.4 for 
different fillings {N — 200 and — 300) and trapping poten- 
tials, z is the coordination number. With vc = 0.2 we mean 
constant trapping Vc/t = const, while vc — mod has parame- 
ters according to Table IVlIIl The dashed line ('[/-0-99(i) ') is 
a power-law fit with exponent —1 within error bars. Compar- 
ison with the homogeneous system at the same temperature 
and on a lattice of the same size is made ('Homog'). Its slope 
is equally -1. 



where n^in and ?imax are the values of the largest and 
smallest value of n{k). 

Deep in the superfluid phase the darkest spot has al- 
most riniin ~ 0, leading to a visibility close to unity, while 
in the Mott insulating phase the contrast between the 
brightest and the darkest spot is almost zero and one 
expects a very low visibility. 

In Fig. [IS] we show the visibilities for a trapped two- 
dimensional Bose-Hubbard model at different densities. 
The calculation was done at a rather high temperature 
(3t = 0.4. Even at this temperature regions with integer 
density and reduced compressibility exist. Precursors of 
the Mott-insulating regions can be seen in the density 
distribution and its variance around U/t ~ 25. Looking 
first at the case of a constant trapping potential, we see 
that for the chosen number of particles a higher density 
leads to a higher visibility at low U/t. The visibility is 
lower at high U/t since the Mott region is larger. Tak- 
ing the steepening of the trap into account we see that 
at low values of U /t the visibility is higher than the one 
obtained with the constant trapping for the considered 
particle number. The reason is the increased number of 
particles with density between one and two which form 
a superfluid edge between an n = 1 and an n = 2 Mott 
region. The n = 2 Mott region is absent in the calcula- 
tion with the constant trapping potential. For the two 
curves labeled 'vc=mod' taking the change in the trap- 
ping potential into account, the visibility is well fitted 
by a logarithmic curve for large values of U/zt > 10 in 
agreement to the finding in Ref. 




FIG. 19: Visibility for a two-dimensional system L = 20 x 20 
as a function of pt at constant particle number = 200 for 
different values of U/t and Vc/t according to Table lyllll 
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this good agreement might be the suppression of the vol- 



FIG. 20: Visibility along the entropic lines of Table |IX] The 
visibility is computed by taking the brightest and darkest 
point of the momentum distribution of the Bose-Hubbard 
model only. 



ume fraction of the superfluid regions, such that only the 
Mott-insulating regions contribute to the decay of the 
visibility. 

In Fig. [19] we show the visibility as a function of 
the inverse temperature for different strengths of the 
optical lattice potential. We see a strong decrease in 
the visibility if the temperature becomes of the order of 
f3t w 0.2 when the quantum regime is left. A comparable 
decrease in the visibility at low temperature by changing 
the interaction strength is only possible when it's very 
large, i.e. U/t = 200 leads to a visibility of the order of 
0.2. However, we should note that at high temperatures 
the single-band approximation of Eq. ([1]) loses its validity. 

We further show the visibility along different isen- 
tropic lines in Fig. [^D] for the same parameters as taken 
in Table IIXI The curves with low entropy (low initial 
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FIG. 21: Visibility y as a function of U / zt along an isentropic 
line 5* = 10 with the same parameters as in Fig. [15] (see Ta- 
ble IVIII|I . We compare the ideal theoretical value obtained 
by taking just the maximum and minimum in the momentum 
profile ('Ipt') with the visibility after averaging over a square 
of length 2n/5 taking the Wannier momentum profiles into 
account ('27r/5 — Wannier'). The error bars are smaller than 
the point size. 

temperature) both have a slope of approximately — 1, 
but the onset value U/t is different. Even for the curve 
with high entropy, corresponding to an initial normal 
state, the data points for large value of U/t seem to be 
consistent with this slope, but differ for a small value of 
U/t (experiments show a constant visibility within error 
bars at low values of U/t). 



B. Comparison to experimental results 

We will now take the trap and the isentropic tempera- 
ture change into account, and compare with experiments. 
This is the hitherto most realistic description done with- 
out numerical approximations, but for the 2D case we fail 
to take the same number of atoms and system sizes as in 
the experiment. Further we ignore time-of- flight collision 
effects. 

We first compare our 2D calculations to the experimen- 
tally obtained visibilities. Up to now we have evaluated 
the visibility using n{k) at specific wave vectors k. In 
experiments instead an average over a square around the 
brightest and darkest spot is taken ( see Ref. [l0| for the 
exact experimental procedure). The size of the area is 
a trade-off between signal and noise. The squares are 
chosen such that the contribution of the envelope of the 
Wannier functions (cf. Eq. (jlSp ) is minimal, but not neg- 
ligible. In Fig. we show the different curves including 
the average and our previous definition of the visibility. 
For low values of the lattice height the curves without the 
Wannier envelope have values very close to one, whereas 
the curve considering the Wannier envelope starts at a 
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FIG. 22: Comparison between the 3d data of Ref. [l3| (blue 
and red; the red data points have approximately 79,000 par- 
ticles and the blue ones 224,000. This would correspond in 
the atomic limit to the presence of a small n = 2 plateau 
for the lower particle number and to a n = 2 plateau with 
almost half the number of particles for the higher number of 
particles[lll].) and the visibility found in our model (labeled 
'Thl' and 'Th2'). The curve 'Thl' corresponds to the same 
parameters as in Fig.[T5]along the adiabatic line S = 10. The 
theoretical data is treated in the same way as explained in 
the text by taking the Wannier functions and averaging over 
a square of size 2tt/5. The curve 'Th2' is a theoretical curve 
computed in the same way as 'Thl' and for the same system 
along an isentropic line S — 5 but for = 50 particles such 
that the maximum occupation in the center of the trap does 
not exceed unity. 



smaller value for the visibility. Fitting the slopes of the 
visibility with a power law gives roughly the same expo- 
nent for the theoretical and the experimental procedure, 
but the quality of the fit is better for the theoretical pro- 
cedure. 

Looking at three-dimensional experimental data, many 
groups find a visibility close to unity until U/zt ^ 6.5, 
where z denotes the coordination number. For larger 
values of U/t the decrease in visibility is well approx- 
imated by - (U/zt)". The Ziirich group Gi, ^ 
finds ly = —1.36(5), while the Mainz group [46[ finds 
ly ^ —0.98(7). The reason for the discrepancy is not fully 
understood. Similarly, the momentum distribution data 
on a lo g-lo g plot of the two-dimensional system studied 
in Ref. |43| were consistent with a behavior of (t/U) and 
thus with a « — 1. 

In Fig. [221 we show a comparison between our theoret- 
ical results for a two-dimensional system and an example 
of the experimental data [46j taken for a three dimen- 
sional system. We rescalc the interaction strength by the 
coordination number taking the mean-field effect between 
the different dimensions into account. We compare two 
theoretical curves with TV = 50 and N = 200 particles. 
In the strong interaction limit this corresponds to a den- 
sity profile which has an n = 1 plateau for = 50 and a 
profile with an n = 1 and an n = 2 plateau for N = 200. 
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The initial temperatures at U/t = 5 are (3t = 1.7 and 
fit = 1.4, respectively. Note that the particle numbers 
and temperature do not directly correspond to the ex- 
perimental ones. For a direct comparison of an exper- 
imental to a theoretical curve, three dimensional simu- 
lation with realistic particle numbers would be needed. 
However, wc find a good agreement in Fig. [21]of the the- 
oretical visibility treated in the way described above with 
the experimental data even up to relatively strong opti- 
cal lattices. This shows that the experimental visibility 
can be explained by just taking the isentropic change of 
the temperature into account. In particular, the system 
does not have to leave the quantum regime to reach the 
low values of the visibility. The drop in the visibility can 
be explained by the formation of a broad Mott-insulating 
region (cf. Fig. [TBI and [T7|) . 

In one-dimensional systems the procedure to get a 
quantity similar to the visibility has to be changed, 
since the interference pattern consists out of stripes. In 
Ref. [11] the superfluid to Mott-insulating transition in 
one-dimensional tubes was considered. The experimen- 
talists extracted the coherent fraction of the atoms by 
taking the ratio of the content in the interference peaks 
and the background. We compared (not shown) the ex- 
perimental data with the theoretical calculations where 
we treated the results for the momentum distribution in 
the same way. We found the same order of magnitude 
and the same qualitative trends ( e.g., the visibility for 
low lattice heights is between 0.5 and 0.6 in both cases). 
A detailed comparison however is hindered by the pres- 
ence of many different parallel tubes in the experiments 
and by the large uncertainties stemming from the diffi- 
culty to fit the very sharp interference peaks. 



VII. SUMMARY 

An interpretation of experimental results on the su- 
perfluid to Mott insulator transition taking finite tem- 
perature effects into account has been given before by 
several authors, with differing conclusions. We have ad- 
dressed the possibility of adiabatic temperature changes 
when ramping up the lattice for a Bose-Hubbard model 
using unbiased and first-principle quantum Monte Carlo 
simulations. We find that the expressions 'heating' and 
'cooling' have to be taken with care, since they strongly 
depend on the unit in which temperature is measured. 

For the homogeneous case, we found some small heat- 
ing in units of t near commensurability in the superfluid 
phase for high temperatures (3t < 0.5. For low tem- 
peratures, there is some small heating for low values of 
U/t < {U/t)c, but the system cools slightly for larger val- 
ues of U jt when we are in the proximity of commensura- 
bility. At very low densities, temperature remains almost 
constant in the superfluid phase. When the density is 
commensurate, the temperature shoots up dramatically 
in units of t when the Mott gap opens in order to keep 
the accessible number of levels constant. This is in agree- 



ment with the findings of previous studies 0, [13, HH . 

In a trapped system however the situation is different 
as already noted in [l^l for hard-core bosons. We find 
that in a one- and a two-dimensional system in a trapping 
potential the main entropy contribution comes from in- 
commensurate regions with low filling n sa 0.5. This is in 
good agreement with the finding for a three-dimensional 
system of the Amherst group where they find that the 
non-commensurate edges between the Mott plateaus be- 
come normal and accommodate the entropy [1^. We 
found that the temperature increased in units of t when 
the size of the commensurate regions is broader than sev- 
eral lattice sites in the case of constant trapping. For 
realistic trapping and densities around one or two, we 
found a temperature increase in units of t, even in the 
weakly interacting regime. However, in units of U tem- 
perature decreases or saturates, and its value lies deep in 
the quantum regime for the commensurate regions. Only 
if the incommensurate regions around filling n sa 0.5 are 
almost totally suppressed the entropy has to be taken by 
excitations in the commensurate regions or regions with 
higher filling, by which the Mott-insulating regions might 
be destroyed. Our result is in agreement with the results 
of F. Gerbier who finds that current experiments easily 
reach the thermal insulator regime (T < T* sa 0.2J7, 
where Mott-like features persist but superfluidity is ab- 
sent), and possibly the quantum region [ist . In contrast 
to the predictions of Ref. [1], we find that for realistic 
paramters no runaway temperature occurs. 

Assuming adiabaticity in current experiments is in 
agreement with theory remaining in the quantum regime. 
In particular, the drop in the visibility of the interference 
pattern can be fully explained within this framework and 
no additional temperature rise has to be taken into ac- 
count. 

Future experiments using the spatially resolved mea- 
surements of the density and higher order correlations 
will be able to confirm the creation of Mott-insulating 
regions as the first evidence was reported in [i^. Unfor- 
tunately, we have seen that the integrated density fluc- 
tuations are not very sensitive to temperature changes 
when one is deep in the quantum regime. 



VIII. OUTLOOK AND CONCLUSIONS 

Our principal assumption that the loading of the lat- 
tice is approximately adiabatic needs to be verified con- 
sidering the dynamical process at finite temperature. It 
will be equally important to extend our investigations 
to different and larger systems. A major goal will be 
to treat the full three-dimensional Bose-Hubbard system 
with a realistic number of particles. Systems with differ- 
ent types of particles will have to be addressed as well, 
since the development of a new energy scale that does not 
scale with JJ will have a negative influence on the pos- 
sibility of reaching the ground state adiabatically. The 
visibility results of theory (49| at low temperature and 
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experiment [451 . |50| are in disagreement for Bosc-Fcrmi 
mixtures. It was believed that temperature is such that 
these systems arc not in their ground state [i^, and a 
study by Cramer et al. [5l| hints at heating effects for a 
weak inter-species couphng. g the ground state. Anti- 
fcrromagnetism and entropy were previously addressed 
for a homogeneous Fermi- Fermi system in Ref [Sll, [s^, HI] . 

This work was motivated by the strongly different 
opinions that existed about temperature effects in the 
Bose-Hubbard model, ranging from a constant tempera- 
ture to a runaway temperature. We have addressed all 
aspects of this discussion in homogeneous and trapped, 
one and two-dimensional Bose-Hubbard systems and the 
results of our work are uni- vocal: compatibility of experi- 
mental with numerical visibility curves (and density pro- 
files) supports that the experimental initial temperature 
of the BEC is deeply in the superfluid phase, and that the 
quantum regime is not left even when adiabatically ramp- 
ing up the optical lattice, even in the presence of a con- 
siderable Mott domain. The temperature scales almost 
linearly with U when the Mott domain is considerable in 
size, but temperature remains a very low fraction of the 
value of U . We have copied the experimental procedure 
of measuring the visibility in our simulations and found 
good agreement. It turns out that the Wannier function 



and the averaging over a small region of the interference 
pattern produce an almost constant shift of the visibil- 
ity as a function of U/t. For deep lattices, the trapping 
potential becomes effectively steeper, and different Mott 
plateaus are forced to form. The main contributions to 
the visibility come from particle-hole excitations giving 
the visibility a slope of —1 as a function of U/t. The 
superfluid volume fraction is too low, but it still absorbs 
most of the entropy, as we could infer from the local den- 
sity approximation. Our results strongly support that 
experimentalists have observed the superfluid and Mott 
phase and their crossover in a trapped system. 
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